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Abstract 

We shall discuss various kinds of geometric bremsstrahlung pro- 
cesses in the spatially flat Robertson- Walker universe. Despite that 
the temperature of the universe is much higher than particle masses 
and the Hubble parameter, the transition probability of these pro- 
cesses do not vanish. It is also pointed out that explicit forms of 
the probability possess a new duality with respect to scale factor of 
background geometry. 



1. Introduction 

Particles in the early universe undergo severe redsliift as a result of the 
cosmic expansion and get decelerated extraordinarily in the comoving coor- 
dinates. Consequently radiation or massless particles may be emitted from 
the decelerated particle. This process induced by the geometry of the uni- 
verse is regarded as a kind of bremsstrahlung. Thus we call it geometric 
bremsstrahlung. This effect may bring about many sorts of decay and emis- 
sion process, which are prohibited kinematically in the flat spacetime. 

Despite of its ease to recognize such a process in the classical mechanical 
meaning, there is a nontrivial aspect on existence of quantum geometric 
bremsstrahlung. Temperature in many interesting situations of the early 
universe is much higher than particle masses and the Hubble parameter. 
Since momentum of the particle is comparable with the temperature, one 
might naively think that we can neglect all the mass parameters even in 
calculation of the transition probabilities. Hence it might be expected that 
the probability of the process is equal to that calculated in massless theories 
in the fiat spacetime, thus exactly zero. However this naive expectation is 
turned out to be wrong by careful analysis. 

In [0 , the first precise analysis has been performed in the four dimensional 
Robert son- Walker universe. They shows that even in a high momentum 
limit there remains the nonvanishing probability of photon emission via the 
geometric bremsstrahlung. 

In this paper we shall give an extended analysis on several kinds of pro- 
cesses of the geometric bremsstrahlung. It includes analyses on the 0^ the- 
ory in arbitrary dimensional spacetimes, the theory with Yukawa interaction 
and the massive vector field theory. It will be shown that the high momen- 
tum limit, or the high temperature limit, does not terminate the geometric 
bremsstrahlung process in such a rather wide class of interactions and in ar- 
bitrary dimension. It is also stressed that a new type of duality can be found 
in the forms of the transition probabilities for renormalizable interactions. 

In section 2 we introduce the spatially flat Robertson- Walker universe 
with past and future asymptotic flat regions and explain why we especially 
consider it. In section 3 we shall give a review on free fields in our model of the 
universe. In section 4 the geometric bremsstrahlung induced by the Yukawa 
interaction is analyzed in detail. In section 5 the (jr" theory in arbitrary 
dimensional spacetimes is surveyed. In section 6 we give explicit forms of 
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transition probability of the process including a massive vector field. 



2. Robertson- Walker Universe 

with Past and Future Minkowskian Regions 

Here let us comment a little bit on our model of the Robertson- Walker 
universe. In this paper only the spatially flat model is argued, but we believe 
that our results will be extended to open and closed models. The geometry of 
the spacetime is expressed by a metric tensor whose form is written as Qfj^^, = 
aitYrj^y, where a{t) is a scale factor function and rj^y is the Minkowskian 
metric. The argument of the scale factor, t, is conformal time and proper 
time T can be defined by dr = a{t)dt. From the metric the Christoffel symbol 
is easily manipulated as follows. 

= 6pd^ In a + 5° 9/3 In a — rip^d"" In a, (1) 

where d" = rj'^^dp. This yields an explicit form of the scalar curvature 
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R=- 



a2 



2(n- 1) — lna + (n- l)(r2-2) I ^Inaj 



(2) 



where n is the dimension of the spacetime. The Hubble parameter is defined 
as usual by the scale factor as 

1 da 1 da 
Hit) = = 

a dr dt 

In later sections we must take account of interaction in the expanding 
universe. Then well-defined asymptotic free field is required to construct S 
matrix elements rigorously. However in general spacetimes the construction 
of asymptotic fields often encounters some difficulties 0. 

For example the universe accompanied with the big-bang possesses a def- 
inite birth time and an initial singularity. Near at the birth time it needs 
more pieces of information (maybe quantum gravity) to specify how to define 
asymptotic in-field. 

The simplest prescription to avoid such problems is to restrict ourselves 
on analysis in the universe equipped with Minkowskian past and future re- 
gions. Clearly this enables us to define both asymptotic in- and out- fields 
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just like those in the flat spacetime. It should be emphasized that phys- 
ically meaningful results derived from this model must be independent of 
artificial detail of the model, that is, how the universe gets into and out of 
the expanding era. It will be found that results obtained later satisfy this 
criterion. 

Owing to appearance of the asymptotic flat regions the scale factor must 
satisfy two constraints. By virtue of the time rescaling invariance one of 
them is expressed, losing no generality, as 

a(t ~ oo) = 1. 

Now let b denote ratio of initial scale factor to final one. Then another 
constraint is written as 

a[t ~ — oo) = b. 

It is worth grasping qualitative behavior of Fourier components of a{t)"', 

/oo 
dta(t)"e*'"*. 
-oo 

Now let ul^^n < < < ■ ■ ■ < c^max denote typical frequencies charac- 
terizing detailed evolution of a(t)" with the following condition satisfied. 

F„(.;f))~0(l). 

(n) 

It is natural to think that all the uul are of order of the Hubble parameter 
in the expansion era. 

For UJ <^ i^min it can be shown that 

F„(a;~0)~-(l-6"). (3) 

UJ 

Note that this relation (j^) holds for rather arbitrary cosmic evolution. 
On the other hand for uj ^ ^^mlx behaves as 

FniuJ ~ oo) ~ 0. 
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3. Free Fields in the Robertson- Walker Universe 



In this section free fields in the universe introduced in section 2 are re- 
viewed for later convenience. 

Let us first review free scalar field in the expanding universe. The action 
with the conformal coupling term is written as 



S. 



scalar 



dJ'xJ\g 



n 



m 



A{n - 1) 



Next let us change the field variable cj) into (p = a(t)"/^ 
is reduced into 



S. 



scalar 



(90)'-i(ma(t))^0 



2 12 



i?j 02 J . (4) 
p . Then the action 

(5) 



This is just a free field action with time dependent mass ma{t) in the fiat 
spacetime. 

Boundary conditions for asymptotic fields can be described in terms of the 
rescaled field 0. It is worthwhile to point out that spatially fiat spacetimes 
possess isometry in the spatial section. Consequently the Fourier transfor- 
mation can be used in equation of motion derived from eqn (|^). Therefore 
what we need is a solution which form is such that 



(6) 



where p is conserved conformal momentum. Its related physical momentum 
is expressed as Pphys = v/ciif)- It is easily shown that this satisfies a 
Schrodinger-type equation, 



m^a{tf 



Uf 



p u 



pi 



(7) 



where p = \'p\- The in- mode function of this equation is specified with the 
following boundary condition, 



1 



~ — oo 



(8) 
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On the other hand, the out-mode function satisfies another boundary condi- 
tion, 



u°/\t ~ oo) 



(27r)"-i2Vp^ + 



(9) 



The in- (out-) mode function may have refiection wave terms in t — » oo(t — >■ 
— oo) induced by the nontrivial potential term, — m^a(t)^. Usually such ex- 
istence of the refiection wave means particle creation from the vacuum state 
in the field theoretical context. Typical energy of the created particle is of 
order of the typical Hubble parameter. Hence the energy can be thought 
much smaller than temperature of the universe. This effect is of some inter- 
est, however it is not our target in this report. In fact we shall concentrate 
our attention on particles with high momentum nearly equal to the temper- 
ature. Consequently the reflection wave term is negligible in the following 
argument. 

Now we can also exhibit results for a soluble example with a step scale 
factor a(t) = bQ{—t) + 6(t). For example analytic form of the in-mode 
function is expressed as follows. 



n^(t < 0) 
v^it > 0) 



1 



(27r)"-i2 + mW 
1 

(27r)"-i2 + m^h^ 



Ae 



where 



A 
B 




(10) 
(11) 



This will be used in section 4. As it should be, the reflection coefficient B 
vanishes in p ^ oo. 

Here we also comment on the WKB amplitudes of eqn(|^). When a high 
momentum condition; ^ m^a^ or p^j^y^ ^ rn^H is satisfied, the WKB 
approximation is validated enough in eqn(0). Then refiection coefficient can 
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be neglected and both in- and out-mode functions are written in the same 
form as follows. 

/{27i)^2E{p,t) 



where E{p, t) ~ ^Jp'^ + w?a{ty. In the early universe situation, this form 
gives us reliable estimation of the wavef unction. 



Next let us review free propagation of spinor field. Writing down the 
action needs viclbcin related with the metric tensor like g^y = ef,e„„. For 
the spatially flat Robertson- Walker universe the vielbein reads 



Then spin connection is obtained from this vielbein such that 



ah 



The action of free spinor field reads 



^spmor 



(i"a;det[e)i] [*i7^V^^ - m^^ 



where 



a I ! 



{7a, 76} = H 



abi 



(7 ah 



d,, - -ujfa. 



ah 



;[7o,76]- 



(13) 
(14) 

(15) 

(16) 
(17) 

(18) 
(19) 



Subsequently by defining a rescale field — a 2 ^ we rewrite the action as 



Sspinor 



(20) 



where we have used the relations; 7*70 — n and 7^7^70 = (2 — 72)7''. Thus, 
like the scalar field, the theory can be also reduced into just a free theory 
with time dependent mass in the flat spacetimes. 
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Next let us introduce mode functions more explicitly for n = 4. Here we 
adopt the standard representation for the gamma matrices, 



1 
-1 



7 = 



a 
-a 



Also we introduce ah^^ and (i^^^ defined by the following equations, 



, ^ ■Pi{h,p) = h\p\i{h,p), 



(21) 



(22) 



where h — ±1 and h/2 is helicity of the particle. Then a^^^ and (5h^p satisfy 
equations such that 



. d . . 

I- malt) 



d"^ d ' 

— + pi + m a{t) + im—a 
aV at 



Oih,p{t) = 0. 



(23) 
(24) 



The in- and out- mode functions for ct/j p- are specified by imposing boundary 
conditions. 



in(±j/, \ 
O^h} '{tr. -OO) 



(27r)32 + mW 



(25) 



- oo) = 4^^±^e-^*V^, 
'"^ ^(27r)32Vp^T^ 



(26) 



where the sign -|-(— ) in a^"^^ corresponds to particle (antiparticle) wavef unc- 
tion. 

For a step evolution like a{t) = bQ{—t) + 0(t), exact analysis is possible 
and for example analytic m(+)-modc function is given by 



(27r)32Vp^T^^ 



where 



2 
1 



jp^^+m?W ^ m(\ — h) 



p2 _|_ ^2 y^p^ + 



(27) 
(28) 



For the spinor field the WKB approximation also can be justified when 
p ^ ma or Pphys ^ fn. Then the following amplitude is obtained. 



(29) 



It will be also useful in the later sections to introduce U and V spinors 
corresponding to particle and antiparticle as follows. 



U{h,p,a) 



^ E{p,a{t)) +ma^ {h,p) 
h^E{p, a{t)) - ma^{h,p) 



, V{h,p,a) 



-h^ E{p, a{t)) — mari{h,p 
^jE{p, a{t)) + mar]{h,p) 



where E{p,a{t)) = ^Jp"^ + rn^aity and rj{h,p) — —ia'^$,*{h,p). Using a polar 
parametrization such that 



p = p{sm 9 cos (p, sin 6' sin 0, cos 9) , 
the explicit forms of ^ and rj are written down as follows. 



e(l,p) = 



?7(1,P) 



e *'^/^cos 



6^-^/2 sin (f) 



COS 



(C 



e(-i,p) 



COS 



(i) 



_g-#/2 
_„i</./2 



COS 

sin (; 



(i) 



(30) 



(31) 



(32) 



Next let us give a review of massive vector field in the four dimension. 
The original equation of motion is written as follows. 

V^F^. + m^A, = 0. (33) 
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where F^y = V^A^ - V^A^ = d^A^, - d^Af, . 



Now using the conformal flatness g^y = a{t)^ri^i,, the equation is rewritten as 

+ m^a^] A^ - d^{dA) = 0. (34) 
The transverse wave solution can be introduced such that 







(35) 



where h = ±1 and n^^^ is a unit vector satisfying p ■ n^^^ — 0. The equation 
of motion requires 



— +p' + m^a{ty 



On the other hand the form of longitudinal wave solution is 

a(^)(t) 



Time dependent factors a^^^ and P^^^ need to satisfy 



p2 _|_ iji^a{ty dt 



Rescaling such that 



jp2 _|_ ^j2q2 

m?a? 



(36) 



(37) 



(38) 



the following Schrodinger type equation also should hold 

2 o 



(f' 2 2 2 Sni^p^ / da^ 

■^+p +ma + ^^2 + ^2^2)2 [-^^ 



p- 



IdM 



p2 _|_ j7|2q2 ^^2 



m = 0. 



In high momentum situation,p » ma, the WKB solution takes the 



following form. 



p 



^{2T^f2pma 

1 » 
— I 
(27r)32p^a 



'i--^ + o(p-^) 

pa dt 



-i J dt 



P+ 



m a 1_ d a 

2p 2pa dt^ 



[l + Oip-'))e 



i J dt 



P+ 



m a 1 d a 

2p 2pa "dt^ 



(39) 
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Thus the longitudinal component is estimated in the high momentum limit 
as 



(27r)32p 
where 



idfj_A — Sofj, exp ^p- x — i J dt ^ + 



m 



1 d^a 



2p 2pa dt"^ 



1 ip-x—i f dt 

A = — e 



2 2 1 j2 
I m a 1_ a a 

P""" 2p 2pa 



ma 

These solutions will be used in section 5 and 6. 

4. Yukawa Interaction in the Expanding Universe 

In this section we shall discuss high energy limit of transition probabilities 
via Yukawa interaction in the four dimensional universe. The interaction is 
expressed by 

Syukawa = ^ j d^X^/^(t){^i<^2 + ^ 2^ l) 

= A y" /x0(§i§2 + ^2^1) (41) 

where \l/i(\l/2) is a spinor field with mass mi{m2) and is a scalar field with 
mass /i. The tilded fields are rescaled ones in section 2. 

This action generates several types of processes. For a particle with 
conformal momentum p and helicity 1/2 to decay into a '^2 particle with q 
and a particle with A;, the amplitude is written down as follows. 

Amp^, = -i\ J d^xe'^^-'^-^'^-^uf'*{t)^<T{hf, g, t)¥;'{h, = 1, p, t) (42) 

where u°^^{t) is an out-mode function of (p and v^™/°"* is an in-/out- mode 
function of ^ and hf/2 is helicity of the created spinor particle. In eqn(|42|), 
we have omitted contribution from Bogoliubov coefficients of vacuum polar- 
ization which is suppressed in the high temperature limit . In the following 
argument we assume that uii < 1712 + yU, so no decay happens in the flat 
spacetime limit. Even in the expanding spacetime the decay is prohibited 
especially in the past and future flat regions. Thus the decay is physically 
interpreted to occur only in the era of the cosmic expansion. 
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We expect that this assumption, mi < m2 + /i, is not essential for our re- 
sults when the decaytime we will obtain is much shorter than that calculated 
ordinarily in the flat spacetimes, E^/ {rnxV fiat)- 

In the expanding universe its scale factor is actually time dependent. 
Hence the meaning of probability per unit time seems ambiguous. There- 
fore wc adopt transition probability itself to get a clear interpretation. The 
transition probability is defined by 

where 

■ I ± r „ ~ in/ out r, ~ ■ / i 

j^m/out ^ dx^ ^ ^O^tn/out^ ^^^^ 

J 

= J^^dxH{uldtU^-dtuluj;). (45) 

Now our purpose is to analyze the high momentum limit {p — ^ oo) of X) I'S'P- 
As a warm up, let us first calculate a simple case with the scale factor 
a{t) = bQ{—t) + Q{t). It is possible to calculate the limit straightforwardly. 
Substituting explicit form of mode functions into the definition of jS"!^, we 
get 



wif'\l/2^hf/2)=lim^\S\ 



p^oo 647r'^ J 



X 



dtr, 2 , 2\/ 2 , TWT^ ^^2~^ 2Ml/4^2(/i/, Q, d — l)Ul{l,p, tt 



rO 



-it(^^yp^+mjb'^-y/q'^+m.lb^-^y{p-q)'^+^l'^b'^^ 



where we have used the fact that the reflection component of the wavefunc- 
tion vanishes in the limit. After some tedious manipulation flnal analytic 
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forms of W turns out to be as follows. 



iy(^;-p)(i/2 ^hf = -1/2) 

wif'\l/2 ^hf ^1/2) 



where 



F{mi, 1712,11) = I dy 
Jo 



327r2 
A2 



l + b\ 1 



In — - 2 



87r2 



1-62 62 

6 In 6^ 



(46) 



1-62 



{l-y){miy + m2) 



IJ?y + mi(l - y) - mfy(l - y) 



F{mi,m2,ijL), (47) 



(48) 



Next we try to obtain W for arbitrary a(i) except a(oo) = 1 and a(— oo) = 
6. Our strategy comes from the fact that approximate energy conservation 
holds for high momentum reactions as seen below. Now we concentrate on 
manipulating contribution from the phase space region where p ^ ma and 
q and k — \p — q\ are of order of p. As a result of this restriction the 
WKB amplitude can be justified for both initial and final mode functions. 
Consequently W is turned out to possess a factor hke 

/oo _ 
dtU2{hf.q.a(t))UAl-lla{t)) (49) 
-oo 



exp 
X 



i J dt (^g2 + rnla{t)^ + ^{p- q)^ + iJ?a{tY - yjp' + mla{ty 



[(p2 + m?a(t)2)(g2 + mla{ty){{p- q)^ + ^Pa{tY)]^'^ 

It is very suggestive for /ly = 1 to roughly estimate it by neglecting masses 
as follows. 

A(/i/ = l) ~ f{p,q) X J dte'^^'i+^P-'^-P'> 

= fip,q)xi27r)'5iq+\p-q\-p). (50) 

Therefore the following relation must be satisfied at least to this leading 
order. 

p ~ g + \p — q\ ^ q -\- \J ip - qY -\- 2pq{l — cos6) (51) 
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where p ■ q = pq cos 9. This is approximate energy conservation law and an 
important clue for us to calculate VF^^. From this "conservation" law only 
phase space region where < q < p and ^ ~ hold can contribute to W^j,^ . 
This fact tempts us to introduce a tiny but not specified constant 9o <^ 1 
and a small constant e satisfying mi, 1712, pe <^ p. Now consider 

only a portion of the phase space with pe < q < p{l — e) and < ^ < ^o- 
Consequently the following expansions are valid; 



Furthermore as a result of the high momentum limit, we can expand 



sin 6 6, 

l-cos^~ V. 

2 




(52) 



(54) 



(53) 



From these useful expressions, we obtain for arbitrary a{t) 



}imY^\S{hf = l)\' = W^,{hf = l) 




Next change the integral variables as follows. 



9 



py, 



mi 

P 



(55) 
(56) 



t 



2p 



(57) 
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Then p dependence appears only in the upper bound of z and the argument 
of the scale factor. In fact the result is expressed such that 



1^^,(1/2^1/2) 



hm — 



1-e 



dy- 



y -Jo 



drjexp 



- [ drj ( 

I Jo \ 



dzz 



(58) 



ml{l-y) ^ 



To evaluate VT^, the following relation is worth proving. 

'2pri 



lim a 



lim a 

p^oo 



2 



be{-r]) + e{7]). 



(59) 



Actually in the limit ^ approaches oo when i] > 0, while — oo when i] < 0. 

From the fact that a(oo) = 1 and a(— oo) = b one can derive easily the eqn 
([591). By substituting ( ^9]) into (|^) we can proceed with calculation of W^*,. 

Because we take p — > oo, Oop/m can be replaced into oo. Hence no 
dependence of 9o remains in the final form of VF^,. Moreover we can take 
e — > because of no appearance of infra-red divergence in the integral. These 
replacements yield 



1^^,(1/2^1/2) 



X 



dy{l-y) 
1 



dzz^ 



1 



z'^ + A{mi; 1712, yU, v) z'^ + b'^A{mi; m2, /i, y) 



where 



A{mi]m2,n,y) 



y' 



-l(\-y)^l^y-y(\-y) 
Till nil 



(60) 



After integrations with respect to y and z, the right hand side ends up to 
the same form of W^^^^p\ 



W^,{hf = l) = wiT'\hf = l). 



(61) 



It is a prominent feature that this relation (|6lD holds for arbitrary a{t). This 
property implies the existence of notable universality of Wis,j^{hf = 1). 
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For hf = —1 the story changes a httle bit. The delta factor in eqn(|^) 
can be estimated roughly as 

A{hf = -1) ~ f'{p,q) X J dta{t)e''^''+^P-'l\-P^ (62) 
~ f'ip,q)F,{AE), (63) 

where AE = q+\p — (l\~p and 

Fi(AE) = J dta{t)e''^'^ . 
This scale factor contribution comes from 

U,{-1, q, a{t))U^{l,p, a{t)) ^ ^ (^I^ + j a{t). (64) 

If < g < p and 6^ ^ 1 do not hold, AE 0{p) ^ ^mlx is inevitably 
satisfied. Therefore, as mentioned in the section 2, Fi ~ holds. Thus 
approximate energy conservation becomes valid again. 

0<q<p, (65) 
^ ~ 0. (66) 

In this case it should be noticed that 6 is not so large that AE becomes 
larger than u^l^. Only phase space region with 



9<9o 



\ 



p 



contributes to VT^i. Estimation of Wx;,-^ for /i/ = —1 is also possible taking 
account of eqn(|67D. Thanks to the high momentum limit the scale factor can 
be replaced the step one as a(t) = bQ{—t) + 6(t). After manipulation we 
can show 

W^,{hf = ~l) = W^f'\hj = -l). 



Consequently eqns (^6|) and (^71) are turned out be correct for arbitrary 
evolution of a{t). Eqn(|59D is the keypoint giving birth of the universality. 
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However in the actual universe p does not take infinite value but finite and 
is of order of the temperature. It should be noticed for large but finite value 
of p that the replacement of a{t) into the step evolution is valid only when a 
high momentum condition like 

2 

^min > (68) 

P 

is satisfied. Here denotes the lowest typical frequency of a{t) and is 
assumed of the same order of minimum value of the Hubble parameter. The 
condition ( |5BD enables us to regard a(t) as 66 (—t) + 0(t). The inequality 
(|68D can be also rewritten as 

1 

— > 
m 

where 1/m represents Compton length and l/uJrnin means maximum radius of 
the Hubble horizon. The factor m/p ~ m/E can be interpreted as Lorentz 
contraction factor. Imagine a particle with Compton length 1/m running 
with high momentum p in the expanding universe. Also suppose that the 
universe begins to expand when the particle reaches a point A and that the 
universe ceases to expand when the particle arrives at a point B. Length 
between A and B can be naively considered of order of the maximum Hub- 
ble horizon (~ l/umin)- Then the particle can get excitation energy from 
the gravitational field only while running between A and B. Meanwhile get- 
ting p larger, the length between A and B becomes Lorentz contracted as 
[m / p]{l / ujrnin) from the particle's viewpoint. Thus when relation in eqn(|69|) 
holds the particle cannot see details of the way how the universe has ex- 
panded. Consequently this yields the above universality. 

The universality is actually very powerful tool to estimate the high mo- 
mentum limit of the transition probability. It can be shown that the uni- 
versality also appear in other kinds of interactions and in any dimension of 
the spacetime. However it should be noted that we must take some care of 
its treatment when the probability grows up to infinity as p oo. Using 
the universality we can naively manipulate the limit explicitly. However, it 
might diverge. In the later section such a phenomenon can be observed ex- 
plicitly and then we should cut p off at a certain large value of order of the 
temperature of the universe. 



m 
p 



(69) 
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It is a rather straightforward apphcation to estimate particle decay 
probabihty into a particle with helicity hi and a \E'2 particle with helicity 
h2- Here we must assume that fi < mi+m2 in order to suppress the transition 
in the flat spacetime. The final forms of the transition probability are listed 
as follows. 



167r2 



1±^1 1 

1-62 "^52 



(70) 



W^{hi = l,h2 = -1) 



87r2 



1-62 



In 6^ G{fj.,mi,m2), 



(71) 



where 



ml{l - - 2mim2y{l - y) + mly^ 

G/i,mi,m2 =/ dy — — ^ . 

Jo mf{l-y)+miy-n'^y{l-y) 



(72) 



Astonishingly it is found that all of these forms of the probabilities possess 
a certain kind of dual symmetry. Interchanging 6 into 1/6 keeps exactly the 
forms in eqns (|i6|), (0), ([70D, (|7lj). This means that physics of the expanding 
universe is related with that of the contracting universe. Note that the 
symmetry is not merely time reversal one because initial one particle decays 
two particles for both cases. This might be related with a kind of hidden 
duality. 



Finally we discuss decay rate of the particles. For example let consider 
decay of ^'i particle. The decay is expected to occur when Wq,-^ ~ 1. For 
small 6, ly^j behaves just as 

A^* A2 1 

^^..^3^1.^. (73) 

where A^* is the number of final modes which contribute to the decay. Thus 
when the universe expands so enough that 

ar 1 le-K^ , ^ 

-L = elFx^ (74) 
6 
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is satisfied, the particle gets decayed. Assuming the radiation dominant 
universe, b can be as follows. 



6 = 



= — 

Consequently we obtain a decay rate; 

1 32tt^ 

Tf^ 2e~l^Hi 



(75) 



(76) 



where = l/2ri is the Hubble parameter at production time of the 
particle. 

5. Decay due to Three Point Vertex 
in Arbitrary Dimension 

The geometric bremsstrahlung process is naturally expected to occur in 
arbitrary dimensional spacetimes. Furthermore the existence is not consid- 
ered to depend on whether couplings of the reaction are dimensionful or not. 
Here we give a rather simple example teaching us this feature. Let us con- 
sider the (f)^ theory with mass m in the n-dimensional spacetimes. The action 
reads 



S 



02 _ IaA'-Z^'I , (77) 



where we have introduced a mass parameter A. Hence the coupling constant 
A is dimensionless. In n < 6 the vertex is merely a renormalizable interaction. 
Meanwhile in n > 6, the interaction is not renormalizable and then A means 
some cutoff scale of the theory. Transition probability of a particle decay 
with conformal momentum p is straightforwardly given in the tree level as 



^ \S\^ = (27r)^("-i)A2A6-'^y r-'q J 



dta 



3-n/2^ * 



It is turned out soon that this expression itself converges for arbitrary fixed 
p. However taking p ^ oo, this grows up to infinity for n > 6 . This comes 
from the fact that the vertex is nonrenormalizable for n > 6. However there 
exists no trouble in the real universe because p does not take infinite value 
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but that of order the universe temperature. Thus in the following argument 
we treat p finite but large compared with m and the Hubble parameter. Thus 
the probabihty is also finite. To get a useful lower bound of the probability 
in p ^ m, we again pick up only contribution from the phase space region 
with Q <9 <9o<^l and pe < q < p{l — e). This implies that 



Y,\S\^>W = (27r)2("-i)A2A' 

rp{l-e) 



6— n 



X 



Jpe Jo 



r(n/2- 1) 

^''"rf^sin"-^^ 



For n < 6 other phase space contributions vanishes in the high momentum 
limit and W gives us the probability itself, not merely a lower bound. 

Here assuming that pe ^ m, we can use the WKB wavefunctions for both 
in- and out- mode functions. Replacing the integral variables like 



Q = 

e = 

t = 



py 

m 

P 
2p 



(78) 
(79) 

(80) 



we obtain 



W 



X 



71 



n/2-1 



1-e 



.n—3 



Ay (27r)"-ir (n/2 - 1 



dy 



l-y Jo 



dzz 



n—3 



By imposing m^/p -C ij^min^'^^ replacement a — bQ{—7]) + Q{7]) is again 
validated. Then the integration with respect to rj can be easily performed 
and we obtain 



W ^ X 



n— 6 



TT 



n/2-1 



X 



Jo 



A J (27r)"-ir (n/2 - 1) A 

1 i,3-n/2 



V) 



dzz 



n—3 



-\ 2 



z^ + C z^ + h^c 



(81) 
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where C = A{m; m,m,y) = {1 — y + y^)/y'^- 

Unfortunately for n = 2 this expression is not suitable and we must 
calculate separately. However the estimation is also possible analytically and 
it results in just 



W{n 



0. 



From n = 3 to n = 6 we can take p — oo and e — in eqn (^). The results 
are as follows. 



W{n = 3) 
W{n = 4) 
W{n = 5) 
W{n = 6) 



Stt V3 
A2 / 4 



- In3 



-)■ 



1 



2Vb 



87r2 
A2 /5 



arctan 



327r2 V4 
A2 



ln3 - 1 



1 

73 
A 

m 



m. 



1 + 



l-fe2 



2v^ 



3847r3 



1 + 6% 1 _ 



Note that all these expressions are invariant under a transformation b —>■ 1/b. 

As mentioned previously W diverges for n > 6 in p — »• cx). Keeping 
a;^""^^-* > m^/p in mind, a lower bound with 6o = m/p is obtained from eqn 
(81). The explicit form with 6 ~ is as follows. 



W{n > 6,6- 0) > 



m 



(27r)"-ir {n/2 - 1) (n - 3)(n - 4)(n - 6) VA6 
Therefore after when the universe expands enough to satisfy 

A 

6 

m 



n— 6 



1/6 - O A 



the particle gets decayed via the geometric bremsstrahlung. 



After all we arrive at a conclusion that in spite of the high tempera- 
ture of the universe the transition probability due to the Yukawa geometric 
bremsstrahlung does not vanish even in the n-dimensional spacetimes except 
n = 2. 
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6. Decay Process including Massive Gauge Field 



In this section we shall survey decay processes including a massive vector 
particle with mass fi. 

Let us first consider that this particle interact with fermions \E'i and ^2 
with mass mi and 1712. The interaction term can be expressed as 

= g j (fx^^{cv + CAl'')l''^2Aa + C.C (82) 

where = e'^Aa = aS^^Aa and Aa = a~^Aa. In the high momentum limit 
we can use again the WKB approximation in the section 2. This enables us 
to calculate W explicitly. 

Let us consider first decay of a transverse component of possessing 
conformal momentum p into ^1 and 'I'2 particles. Just as the Yukawa inter- 
action case, we only take account of contribution from restricted phase space 
region as follows. 

Wt (h, h 



hm / dk- — - / dee 

p^oo 2>2-n-^ Jpe p{p — k) Jo 



X 



dtUi{h,p- k, a)(cy + CAj^)j''V2{h, k, a)e^p{p, a) 



U "*(^ 2^ 2(F^fc) 2k- J 



where 



eW(p,a) = ±-^[0,l,±i,0]. (83) 

For any combination of {h,h), it can be straightforwardly shown that Wt 
converges in p 00 and e — 0. After a straightforward manipulation we 
obtain the following results. 
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where 



H 



1 \cv + CArmjy^ + \cv - CArm|(l - yf + 2(|cyp - |cAr)mim2Z/(l - y) 



Jo 



mly + mi(l - y) - iJ?y{l - y) 



Wr(/i = = -1) = 0. 



(84) 



Wt(/i = 1, /i = -1) 



247r2 



1 + 6% 1 



iyT(/i= 1) 



247r2 



1 + 62 1 ■ 

r^^"6^-' 



To get W for longitudinal component decay is achieved by taking the 
photon helicity vector as 



k 



,.^,._li,o,o,-ii. 



(85) 



After some manipulation Wl is in fact converged with h — 1 and h — —1 as 



WL{h^l~h^-l) = 



87r2 



1 + 



61n62 
1^ 



where 



K{mi,m2,iJ,,cv,CA) 



X 



mly + ml{l - y) - ii'^y{l - y) 

cvirrii - 1712) - CA(mi + 777,2) / y 
mi ' 



i-y 



Cy(777i - 7772) + 0^(7771 + 7712) /l-?/ 
77721 



y 



-2{cv + CA)n^y{l-y) 



(86) 
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WLih = -l,h = l) 



87r2 



1 + 



6 In 6^ 



1-62 



K{mi,m2,n,cv,-CA). (87) 



While W{h = ±1, h = ±1) is given as 



Wiih = ±l,h = ±1) 



g'^ \cv (mi - 1712) ± CAimi + 1712] 



167r2 



1-62 ^52 



It is also possible to estimate W when a \E'i particle with momentum p 
and helicity 1/2 decays into a \E'2 particle and a particle. Substituting 
the wavefunction forms in the section 2 into W it is written down explicitly 
for emission of transverse component such that 

W^, {hi;hf,T) 

= lim ^— dk—^^—- I dee 



327r2 Jpe p{p - A;) Jo 



where 

~ ~2k~ ^ 2{p - k) 2^' 

For cases with hi = 1 and hf = 1 our WKB results include an infra 
divergence. 

,y.(l;l,l)^|l|ev-c.p(i±|.ni-2)/|^. (88) 



ly* (1; 1, -1) = £;\cy - ([^ - 2) /' 1^(1 - yf. (89) 

In spite of the appearance in eqns (pSf ) and (^) of infra-red divergence 
with respect to pe, we believe that exact form of the probability converges as a 
result of the existence of natural infra-red cut off fi. This cut off works only in 
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the low momentum region where the WKB approximation is invahd. Because 
we have discussed only where the WKB approximation is valid, our 

treatment is supposed to include such a superficial infra-red divergence. From 
this point of view, it can be naively expected that e ~ ^/p. 

However if /i is exactly zero, no natural infra-red cut off appear. This 
infra-red divergence really exists owing to the existence of the soft particle. 
The number of soft particles is counted only to the accuracy of observations. 
Thus for very soft particles one cannot discriminate between virtual emission 
which contributes particle mass operator and real emission. This enables us 
to add a part of the mass operator contribution to the emission probability. 
Then the theory is hoped to have no real infra-red divergence in physical 
interpretation, similarly to the flat spacetime case|^. 

For other cases of transverse emission, W is converged and flnal results 
are given by 

,2 



R{mi,m2, fi), (90) 



47r2 
where 

'•I y{\cv - CAl^ml - 2(|cyp - \cA\^)mim2{l - y) + |cy + CApm?(l - yf) 



R = I dy 
Jo 



mly + iJ?{l-y)-miy{l-y) 

iy^,(l;-l,-l) = 0. 



For the emission of longitudinal component we get 



lim / dk— -, , 

327r^ Jpe p{p - k) Jo 



dee 



X 



I dtU2ihf,p- k, a){cv + CAl'')l^Ui{KJ,a)e*}''\k,a)e'^^'^^''^^^^^'' 



Note that ^^^.^(1 : possesses an infra-red divergence with respect to 

e ^ 0. 



-2 /. 61n62 



9 
87r2 



1-62 



(5(mi,m2,/i), 
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where 



Q 



X 



dy y^^~y^ 

-y) + mly - mfy(l - y) 
cy (mi - 777,2) + (mi + 777-2) "^^2 



cv(mi - m2) - CA(mi + m2) 



/i 



-mi 



-2{cv - ca)ij> 



(91) 



and it is expected that e fi/p. 

Finally it can be shown that iy^(l;— 1,L) converges and takes the fol- 
lowing form. 



W^{1;-1,L) 



327r2 



cv{mi — m2) - CA(mi + m2) 



1 



52 1 
In — 



1-62 62 



Let us next comment on a three point interacting between the massive 
vector field with mass /i and complex scalar fields </?i and </?2 with mass mi 
and m2 . The interaction is described by 

S^^A = 9 j d^x^f^Aiji {(piV(p2 - V(pi(p2) + c.c. 

We can also give in this case explicit forms of W in each process as follows. 
For decay of the transverse component , 



^(T) ^ _9_ 
487r2 

For the longitudinal component, 
^2 61n62\ 



1 + 1 
1-62 "{,2 



9 



87r2/i2 

r.1/2 



1-62 



(mi - m2)^ + 4:fx^{mf - ms) + A/x'^x'^ 



+ 2m2 — /x2 + 4(mi — m^x + A/j.'^x^ 
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For transverse vector particle emission from a scalar particle, 



1 + 52 1 



87r2 6 

For longitudinal vector particle emission out of a scalar particle, 

6 In 6^ \ /■! dy (1 — y) [{ml — m\)y + (2 — |/)/x^]^ 




1-b'^JJs y ii^{l-y) + mly-m\y{l-y) 
6 In 6^ \ 



.-6V 

Thus probability of geometric bremsstrahlung does not vanish, even in- 
cluding the vector field. 
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